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^ ! Abstract 

We explore the discovery prospects for B'B' pair production followed by direct decays B' — )• hX, 
where B' is a new quark and X is a long-lived neutral particle. We develop optimized cuts in 
the {niBi ,mx) plane and show that the 7 TeV LHC with an integrated luminosity of 1 (10) fb~^ 
k> \ may exclude masses up to m^/ ~ 620 (800) GeV, completely covering the mass range allowed 

H \ for new quarks that get mass from electroweak symmetry breaking. This analysis is applicable to 

other models with hhlJ^T signals, including supersymmetric models with bottom squarks decaying 
directly to neutralinos, and models with exotic quarks decaying directly to GeV-scale dark matter. 
To accommodate these and other interpretations, we also present model- independent results for 
the hblf^T cross section required for exclusion and discovery. 

PACS numbers: 14.65.Jk, 13.85.Rm, 95.35. +d 



I. INTRODUCTION 

This is an exciting time for TeV-scale colliders, with experiments at the Tevatron and 
Large Hadron Collider (LHC) collecting data at unprecedented luminosities and energies. In 
this study, we explore the prospects for discovering new physics through B' B' production, 
followed by the direct decays B' — )■ 6X, where B' is a new down- type quark (with electric 
charge qs' = —\) and X is a long-lived neutral particle, leading to the signal bblpT- This 
study complements our previous study of T'T' production (where T' is an up-type quark 
with qt' = |), followed by the direct decay T' — )• tX, leading to the signal tilpT pj- 

The possibility of new physics leading to heavy flavor signals is, of course, well- 
appreciated, but such signals are usually accompanied by other visible particles from multi- 
step cascade decays. The direct decays considered here are much less studied, but are 
well-motivated from many perspectives. The gauge hierarchy problem, for example, moti- 
vates top and bottom partners, new particles that cancel the radiative contributions from 
bottom and top quark loops to the Higgs boson mass. The fine-tuning in such models is 
generally reduced when these partners are light, making it natural that such particles are 
among the lightest new particles and decay without cascades. The canonical example is su- 
persymmetric models with top and bottom squarks that are lighter than the other squarks 
and decay directly through h — )■ 6x5 ^^^ t — )■ txi, where Xi is the hghtest neutralino. 

Dark matter provides another general motivation for the signals we consider. In many 
models of dark matter, the dark matter particle X is the lightest particle charged under 
an exact symmetry, "dark charge," and it may scatter off normal matter through processes 
Xq — 7- Q' — )■ Xg, where Q' is another new particle. This possibility is especially motivated 
at present by the possibility that such signals may in fact have been seen at DAMA |2|, 
CoGeNT |3|, and CRESST |^. In such scenarios, the Q' particles are necessarily colored 
and have dark charge; they can be produced through qq/gg — )■ Q'Q' and decay directly 
through Q' — ;■ qX. Although these decays may be to any quark flavor, decays to b and 
t are realized in concrete scenarios with WIMP dark matter, WIMPless dark matter, and 
asymmetric dark matter, as we review in Sec. [TTl Such models are also, of course, much 
more amenable to study at hadron colliders than those in which the decays are solely to 
light quarks. 

In a previous study [l[, we investigated the collider reach for up- type quark pair produc- 
tion qq/gg -^ T'T', followed by T' — )■ tX. In this work, we analyze the pair production of 
down-type quark B'B' with subsequent decay of B' — )■ bX at both the Tevatron and 7 TeV 
LHC using the Madgraph/MadEvent/Pythia/PGS4 packages. For relatively small values 
of mx, we will find that B' masses up to 440, 460 and 480 GeV may be excluded given 
integrated luminosities of 5, 10, and 20 fb~^ at the Tevatron, respectively. This reach is 
greatly enhanced at the 7 TeV LHC: with an integrated luminosity of only 100 pb~^, the 
7 TeV LHC's 95% CL exclusion reach is comparable to that of the Tevatron with 20 fb~^. 
The whole region of rriB' allowed by Yukawa coupling perturbativity can be explored with 
1 fb~^ of data, and, with 10 fb~^ of data, B' masses up to 800 GeV may be excluded. The 
7 TeV LHC also has great potential in terms of B' discovery: 3a discovery contours reach 
B' masses of 540 and 700 GeV for integrated luminosities of 1 and 10 fb~^, respectively. 

We also present model-independent results for collider reaches as a function of the bXbX 
production cross-section, with nix = 1 GeV. From these, for any theoretical prediction 
for a{B'B') x B{B' — )■ 6X)^ as a function of rriB', one can easily determine the expected 
exclusion and discovery reaches in ttib' ■ Our results may therefore be applied to other models 



that give rise to the bb^T signal. 

In Sec. ini we discuss models that yield the bbfiT signal and their implications for the B' 
and X masses. In Sec. IIIII we discuss existing bounds on these scenarios. Our simulation is 
described in Sec. HVl The results are presented in Sec. |V]and summarized in Sec. |Vll In the 
Appendix, we list cross sections at the Tevatron and 7 TeV LHC for standard model (SM) 
backgrounds and several benchmark points after various levels of cuts. 



II. MODELS 

We now discuss models that yield the bbJ^T signal and their implications for the B' and 
X masses. We begin with the familiar examples of supersymmetry and universal extra 
dimensions (UED), where the spins of the B' and X particles are {Sb',Sx) = (0, |) and 
(i, 1), respectively. These are model frameworks in which existing searches have been carried 
out and the X particle is WIMP dark matter. We then discuss models with WIMPless and 
asymmetric dark matter, where the spins are {Sb', Sx) = (^, 0). In these models, X is again 
dark matter, but the light mass range mx ~ 1—10 GeV is particularly motivated by currently 
claimed signals. Strictly speaking, our analysis is valid only for the spin assignment of the 
WIMPless and asymmetric dark matter cases, but as described below, it is also applicable 
to the other scenarios with minor modifications. 



A. Supersymmetry with Light Bottom Squarks 

Supersymmetric models yield the bb^r signal when bottom squark pair production is 
followed by direct decays b — )■ bxi, where Xi is the lightest neutralino, an excellent dark 
matter candidate p, ISJ . Squarks are often assumed to decay through cascade decay chains. 
In contrast to other squarks, however, bottom (and top) squarks have their masses reduced 
by the impact of large Yukawa couplings on renormalization group evolution, and their 
masses are also directly constrained by naturalness. There are therefore reasons to expect 
the bottom squarks to be relatively light and decay directly to the lightest supersymmetric 
particle, even if other squarks are heavy and decay through cascades. 

The salient features for this analysis are 

• The signal arises from bottom squark pair production followed by 6 — > 6xi- The 
decaying particle is a scalar, in contrast to all other examples discussed below. 

• The b mass is only constrained by direct searches discussed in Sec. IIIH which require 
mi ^ 230 GeV for small ttl^o. The mass limit on m^ becomes much weaker for small 
mass splitting m^ — m^o. 

• The neutralino mass satisfies m^o ^ 47 GeV, assuming gaugino mass unification (7[. 
Without gaugino mass unification, there is no lower bound on the neutralino mass, 
which in general may be anywhere in the range ^ m^o < mr. 



B. Universal Extra Dimensions 

UED models give the desired signature, where the new down-type quark B' is identified 
with the Kaluza-Klein bottom quark b^. The b^ can be pair-produced and then decay directly 
to Kaluza-Klein hypercharge gauge bosons B^, which may be WIMP dark matter [8|, |9|. 



The salient features of this model are 

• The signal arises from b^ pair production followed by b^ — )■ bB^. 

• The b^ mass is set by the size of the extra dimension. Its mass is constrained only by 
the direct searches discussed in Sec. IIIIl which require rrihi ^ 440 GeV for small tti^i . 
If mfii ~ rrihi , the mass limit is much weaker. 

• The size of the extra dimension typically sets the size of all the Kaluza-Klein particles, 
and so UED spectra are typically expected to be compressed relative to supersym- 
metry. One therefore expects tti^i ~ Tnf,i, but B^ masses anywhere in the range 
^ m^i < rrihi are experimentally viable. 

C. WIMPless Dark Matter 



In WIMPless scenarios [10[, dark matter is in a hidden sector. These scenarios have the 
notable feature that the dark matter candidate automatically has approximately the correct 
relic density, regardless of the candidate particle's mass. 

The WIMPless dark matter particle X couples SM quarks to new quarks through Yukawa 
interactions 

V = XlXQ'^qLi + X1Xf'j,um + XfXB'j^dm , (1) 

where X is assumed here to be a complex scalar charged under a discrete symmetry, q'[^ = 
{uLi, dLi), Ufii, and dju are the SM quarks of generation i, and Q'^^ = (T[, 5^), Tjj, and B'p^ 
are the new quarks, also charged under the same discrete symmetry as X. 

In general, the Yukawa couplings can couple X to any of the SM generations, subject 
to flavor constraints [ll|. Although it is difficult to know what a "natural" flavor structure 
for new quark couplings should be, one reasonable possibility is that new quark couplings 
follow the observed Yukawa hierarchy and couple new quarks dominantly to third generation 
quarks. In fact, for C(l) Yukawa couplings, WIMPless models with dark matter coupled 
to 3rd generation quarks may explain the reported dark matter signals from DAMA and 



CoGeNT 12 



For the purpose of the analysis presented here, the salient features are 

• The signal arises from B' pair production followed by the decay B' — )■ bX. Dark charge 
conservation forbids the cascade decays B' — )■ Wq and B' -^ Zq, and the possibihties 
B' —7- dX, sX are excluded by hand. 

• T' and B' are new quarks that get mass through electroweak symmetry breaking, with 
fnT',B' = Xt',b'v/v2, where v ~ 246 GeV. Yukawa coupling perturbativity requires 
X%, ^ 47r, which implies the upper bound m^/ ^ 600 GeV. For small mx, direct 
searches IsHlSJ place a lower bound, m^/ ^ 440 GeV. This lower bound is weakened 



for larger dark matter mass. 
• For X to freeze out with the correct relic density and preserve a key motivation for 
WIMPless scenarios, it cannot be extremely light, but the range 10 MeV ^ mx < itib' 
is allowed. However, if one hopes to explain the DAMA [2| and CoGeNT |3| anomalies 
and be marginally consistent with stringent exclusion bounds from CDMS [l6| and 
XENONlO/100 [l3l, hght masses with mx ~ 7 GeV are preferred. 

There are several other models that share the basic features described above. One well- 
known example is little Higgs models, where the new quarks are not 4th generation quarks. 



but instead arise from the extra degrees of freedom needed when the gauge symmetry SU(2)i: 
is enlarged. Unhke the WIMPless case, the mass of the new quarks is not generated by 
Yukawa couphngs to the SM Higgs, and thus is not bounded from above by perturbativity. 
These quarks may also be charged under T-parity, decaying to SM quarks plus dark matter 
(the lightest particle charged under T-parity) jl8| . Another example is provided by a recent 
set of models in which U(1)b is a gauge symmetry 19|. In this case, the new quarks are 



4th generation (possibly mirror) quarks that are added to cancel the U{1)b mixed anomaly. 
They are charged under a new U(l) global symmetry, and the dark matter is a scalar which 
is the lightest particle charged under this global U(l). 

D. Asymmetric Dark Matter 

Another class of models predicting the bbJ^x signature are models of asymmetric dark 
matter arising from hidden sector baryogenesis |20|. In this framework, sphalerons of a 
hidden sector gauge group generate a baryon asymmetry by producing exotic quarks J2l|. 
The B' is the down-type exotic quark, which decays to dark matter and right-handed bottom 
quarks, B' — )■ bj^X. Our analysis is directly applicable when the dark matter particle is a 
scalar.^ In this model, the number density of the dark matter candidate is determined by the 
baryon number density. Moreover, the dark matter multiplet and 6-quark multiplet are both 
chiral under a 11(1)^3/? gauge group; both mx and nib are determined by the symmetry- 
breaking scale of U(l)'jn3_R, so one expects mx ~ rrii,. As a result, this asymmetric dark 
matter candidate naturally has approximately the correct relic density. 

For the purposes of this analysis, the salient features of the asymmetric dark matter 
model are 

• The exotic down-type quark B' is not a 4th generation quark, and its mass is not 
generated by electroweak symmetry breaking. As a result, there is no upper bound on 
its mass, which is only constrained by direct searches to satisfy m^' ~ 440 GeV for 
small mx- 

• The asymmetric dark matter should have mx ~ 1 — 10 GeV to correctly explain the 
relic density. But a dark matter mass m,x ~ 7 GeV is preferred to explain reported 
direct detection signals. 

• The exotic quark need not be down-type, but if it is, it necessarily decays through 
B' -^ bX. 



III. CURRENT COLLIDER LIMITS 

In this section, we summarize existing constraints on the bbJ^x signature and related 
channels. All bounds quoted below are 95% CL constraints. 

As noted in Sec. HI we are interested in B'B' production followed by direct decays B' — >■ 
bX. These differ from searches for conventional fourth generation quarks, which we denote 
by the lowercase b' and t', which typically decay through cascades. Nevertheless, we begin 
with these as a useful reference point. Searches for f and b' has been performed at both 



^ In general, the dark matter candidate can also be spin-i for this model. In this case, the relevant signal 
is pair production of the down- type exotic squark, B' , followed by the decay B' -^ bjiX. 



Run II of the Tevatron 122-25 and the LHC 26, 27 . These searches assume that the fourth 



generation quarks couple to the first three generations and decay through 5', t' — )■ qW . For 
6', the most stringent result at present is from CDF searches for h'h' production followed by 
y — )■ tW . The lack of an excess in 4.8 fb~^ of data implies rrift/ > 372 GeV [23|. For t', 
CDF finds no signal for t'i' production followed by t' — )■ hW in 5.6 fb~^ of data, implying 
rriti > 358 GeV \24\, and a D0 search for t'i' followed by t' — )■ qW in 4.3 fb~^ of data requires 
mt> > 296 GeV [25[. At the LHC, null results from CMS searches for b' -)■ tW using 34 pb"^ 
of data imply mb' > 361 GeV [26|. ATLAS analyses of t' or b' — )■ qW in 37 pb~^ of data 
imply rrit'fii > 270 GeV [27]. We stress again, however, that the limits of this paragraph do 
not apply to the B' and T' searches we consider here, as decays B' — > qW and T' — )■ qW 
are excluded by dark charge conservation. 

As noted in Sec. Ill Al however, the bbJ^T signal is produced in the case of supersymmetry 
with bottom squark pair production followed by 6 — )■ bxi, where Xi is the lightest neutralino. 
Both CDF and D0 have searched for this signal. D0 finds no excess in 5.2 fb~^ of data, 
requiring mr > 247 GeV for m^o = and excluding 160 GeV < rrii < 200 GeV for m^o = 
110 GeV 131]. The corresponding CDF result for m^o = using 2.65 fb~^ of data is mr > 



230 GeV IJ]. Taking into account only the difference in B'B' and bb* cross sections, the 
D0 bound m^ > 247 GeV implies rriB' ^ 365 GeV. As we will see in Sec. El an optimized 
collider analysis would imply m^/ ^ 440 GeV when cuts similar to those of the D0 sbottom 
searches are applied to B'B' pair production. Signal significance would be reduced by a 
trials factor associated with the choice of optimum cuts, however ._ 

A search for gluino pair production with g ^ bb followed by 6 — ?■ 6xi has been carried 



out at CDF using an integrated luminosity of 2.5 fb ^ [15|. Candidate events were selected 



requiring two or more jets, large fix, and at least two 6-tags. Using neural net analyses, CDF 
found rrig > 350 GeV for large mass splittings Am = nig— mi > 80 GeV, and about 340 GeV 
for small Am ~ 20 GeV. The result for the case of small Am, where two 6-jets are soft and 
sometimes missed, can be applied to the B'B' search and imply roughly m^/ > 370 GeV. 

SUSY searches by the ATLAS Collaboration at the LHC with y/s = 7 TeV and 35 pb~^ 
luminosity (28[ studied the process of gluino and sbottom pair production with g — ;■ bbi and 
bi — )■ bxi- Events are selected by requiring large J^t and at least three jets, of which at least 
one is 6-tagged. For m^^ < 500 GeV, this search implies m^ > 590 GeV. This limit also 
bounds B'B' — )■ bbjlpT-, where the additional jet results from QCD radiation. It is, however, 
not straightforward to obtain the mass limit on mB' without detailed collider analyses. 

Although we focus here on B' production, if the B' mass is generated by electroweak 
symmetry breaking then the mass difference between B' and T' is constrained by electroweak 
precision data to be less than about 50 GeV (29|. Therefore, bounds on T' production are 
also relevant. The discovery prospects for T'T' production followed by direct decays were 
evaluated in Ref. [l|, and we summarize current bounds here. 

A CDF search for T'T' — )■ tilf^T in the semi-leptonic channel in 4.8 fb~^ of data implies 
m,Ti > 360 GeV for m^x < 100 GeV [30|. SUSY searches for stop pair production at the 
Tevatron followed by t — )■ biu also imply bounds on T'T' — )■ ti^T when both tops decay 
leptonically. Null results from searches at both CDF 3l| and D0 |32] in 1 fb~^ of data 
imply TTij- > 180 GeV for m^c ^100 GeV. Accounting for the difference in tt* and T'T' cross 
sections, this implies the bound m^' ^ 263 GeV. The CDF Collaboration has also reported 
a search for top squark pair production based on an integrated luminosity of 2.7 fb~^, using 
the purely leptonic final states from pp — )■ iiil, followed by ti — )■ bxt ~^ bii^x^ |33|. The 
data are consistent with the SM background, leading to the constraint r/ij^ ^ 150 — 185 GeV, 



where the exact hmit depends on m^o, 'rn^±, and B{xf -^ l^uxi)- 

Stop pair production (either direct or v ia g luino decay g ^ tti) has been searched for 
at ATLAS with ^i = 7 TeV and 35 pb'^ 0. Assuming B{ii -> hxf) = 1 and B{xf ^ 
XiW^*^) = 1, searches have been performed in the semi-leptonic channel with 1 lepton, 2 
jets (including one 6-jet) and large ^t- For 130 GeV < m^ < 300 GeV, this search implies 
rrig > 520 GeV. Cross sections for gg + titl around 8 to 40 pb have been excluded for 
400 GeV < rrig < 600 GeV. It is less straightforward to translate this search limit to the T' 
case, given the very different cut efficiencies of the dominant gg process. 



Finally, a recent search at ATLAS with ^/s = 7 TeV and 35 pb~^ [3J| for pair production 
of fermionic top partners decaying to a top quark and a long-lived neutral particle gives 
a mass limit of rrix' > 275 (300) GeV for mx < 50 (10) GeV. This limit can be directly 
applied to the case of T' pair production followed by direct decays. 



IV. COLLIDER ANALYSIS 

A. Signal and background simulation 



Both signal and backgrounds were simulated using MadGraph/MadEvent 4 [35[ and 
passed through Pythia 6.4 [36| (with p^-ordered showers) for parton showering and 
hadronization. We used the CTEQ6L1 parton distribution functions [STJ and the factor- 
ization and renormalization scales were set to rn^ + p"^ of the massive particles produced. 
Detectors were simulated with PGS4 |38| using the Tevatron and ATLAS detector cards for 
Tevatron and LHC simulations, respectively, as provided by MadGraph/MadEvent. 

Our signal process is pp/pp — ?• B'B' + jets with B' — )■ bX, using matrix elements with jet 
matching for up to 2 jets and the decay at matrix element level, giving a signal of two 6-jets 
and missing energy, plus possible associated jets from QCD radiation. We generated events 
at grid points in the {mB',mx) plane with 25,000 events per grid point. 

The main backgrounds to this process are Vr^+ jets, Z+ jets, and ti production. The 
former two were simulated with up to 3 jets coming from the matrix element and the latter 
with up to 1 jet from the matrix element, to ensure that the backgrounds were properly 
modeled with respect to the jet cuts used in the analyses. Also the single top background 
and diboson background were simulated and found to be negligible, as expected. The back- 
grounds have been compared to similar Tevatron |39| and ATLAS 40|-|42| analyses, with 



agreement at the expected 20% level. The exception is 6-tagging, where the experimental 
efficiency for the Tevatron D0 detector is better than that given by PGS. We have therefore 
applied a correction factor to our Tevatron 6-tagging efficiency to reproduce the efficiencies 
quoted in the experimental analyses. 

The choice of cuts to distinguish signal from background is guided by a few key features. 
SM backgrounds exhibit ^t either because of neutrino production or jet energy mismea- 
surement. The first source is suppressed by a lepton veto, which rejects processes involving 
W — )■ iu. The second source can be effectively suppressed by a combination of minimum ^t 
cuts (which suppress multi-jet backgrounds) and the requirement that J^t not be aligned 
with any energetic jets. The alignment cut is imposed in terms of a minimum angle between 
the missing transverse energy and any of the selected jets, A0inin(^T, jets). Since ^t can 
also be mismeasured in events with significant transverse momentum in low-energy jets and 
leptons, such events can be removed by cuts on the quantity A = {$t — ^t)/{,$t + -^t), 



where ^t = \ Z^jets Pr I is the magnitude of the vector sum of the transverse energy of all 
jets. 

Moreover, signal events result in the pair production of two massive objects, whose decays 
to 6-jets and invisible particles are roughly uncorrelated. These events are thus expected 
to have several objects with large pt, with most of the transverse momentum carried by 
the two leading jets, which are not expected to be back-to-back. The J^t is also expected 
to be comparable to the pt of the other objects. For the Tevatron, the presence of objects 
with large pt is measured by the kinematic variable Ht = SjetsPr' while at the LHC it is 
measured by the kinematic variable M^a = $t + Y.j^...jiVT- 

At the Tevatron, the requirement that the leading jets not be back-to-back is measured by 
the kinematic variable Uj^j^i defined as the angle between the two leading jets in the trans- 
verse plane. At the LHC this requirement is imposed in terms of the transverse sphericity, 
5*7-. If Ai,2 are the eigenvalues of the 2x2 sphericity tensor Sij = Y^kPkiV^^ for all selected 
jets, one defines St = 2X2 1 {Xi + A2). QCD backgrounds are dominated by back-to-back jet 
configurations, for which St ~ 0. 

At the Tevatron, the requirement that most transverse momenta be carried by the two 
leading jets is measured by the kinematic variable Xjj = {pj^ + pj:)/ Ht- At the LHC, 
the requirement that the missing transverse energy be comparable to the momenta of other 
objects is measured in terms of / = -^r/^eff- 

Note that we have chosen not to take into account next-to-leading order K factors in our 
analysis. We consider this conservative, in the sense that next-to-leading calculations tend 
to increase signal and background (both QCD and vector boson + jets) with a similar factor 
(~ 1.2 — 1.5), and therefore tend to increase the significance of the result and improve the 
exclusion and discovery regions |43|, |4^ . 



B. Tevatron Cuts 

For Tevatron searches, we begin by imposing precuts that are similar to the cuts required 



in Ref. |39[ for sbottom searches. With these precuts, we require: 



lepton with \r]e,\ < 2.0 and p^^^ > 15 GeV and pJp > 10 GeV. 
2 or 3 jets with \rij\ < 2.5 and p^ > 20 GeV. 



C(jij2 — 165°. 



^T > 40 GeV, ^t/GcV > 80 - 40 x A0„,in(|^T, jets). 

At least two jets, including the leading jet, are tagged as 6-jets. 

A0inin(-^T, jets) > 0.6 rad. 

-0.1 <^<0.2. 

Xjj > 0.75. 

In addition, for each grid point in (m^/, nix) space, we consider the following cuts and 
choose the combination that optimizes the signal's significance: 

• pi^ > 50, 80, 100, 150 GeV. 

• ^T> 100, 150, 200, 250 GeV. 

• Xjj > 0.9. 

• Ht> 150, 220, 300 GeV. 



Note that some of the final cut combinations are redundant (e.g., Ht < 'Pt^'Pt)i resulting 
in a total of 160 naively non-redundant combinations. 

Comparing our cut efficiencies for sbottom pair production with those listed in Ref. [39| 
for two signal benchmark points (m^ , m^o) = (240 GeV, GeV) and (130 GeV, 85 GeV), we 
found reasonably good agreement, except for the 6-tagging efficiency for two jets, which was 
underestimated by about a factor of 2 in PGS4. We therefore increase our cut efficiencies 
by a factor of two (for both the signal process and the ti background) to account for this 
underestimation caused by our using PGS4 for detector simulation. 

C. LHC Cuts 

For the LHC, we adopt the following precuts based on the cuts designed for inclusive 



SUSY searches for lepton, 2-3 jets (including 1-2 6-jets) at the LHC |40H42l: 

lepton with |?7e,^,T| < 2.5 and p^^'^ > 20 GeV. 

2 or 3 jets with |r/j| < 2.5 and pj]^ > 100 GeV, pf' > 40 GeV, p^^^*"^ = 30 GeV. 

^T > 80 GeV. 

/ = ^r/Mes, f > 0.3 (0.25) for 2-jet (3-jet) events. 

A0mm(-^T, jets) > 0.2 rad for all selected jets. 

Transverse sphericity St > 0.2. 

At least one selected jet is tagged as a 6-jet. 

The cuts on the transverse momentum of jets are chosen to satisfy the trigger require- 
ments, as well as to reject a sufficient amount of QCD jet background. 

In addition, for each grid point in {rriB', mx) space, we consider the following cuts and 
choose the combination that optimizes the signal's significance: 

• p^ > 150, 200, 250, 300 GeV. 

• $T> 100, 150, 200, 250, 300 GeV. 

• Meff > 250, 300, 400, 500, 600, 700 GeV. 

Note that some of these final cut combinations are redundant (e.g., when Meg < p^ + 
Pt + -^t); resulting in a total of 104 non-redundant combinations. 

V. RESULTS AND DISCUSSION 

Figure [1] shows the IpT and Ht distributions for three benchmark points {mB'-,mx) = 
(200, 1), (300, 1), (400, 1) GeV as well as the dominant SM backgrounds at the Tevatron 
after precuts. W{£h')jj becomes the dominant background after precuts, and Z{uu)ii and 
semileptonic ti background are relatively large as well. While the differential cross section 
distributions for the SM backgrounds drop quickly with increasing ]pT and Pt-, the dis- 
tributions for the signal typically extends to much larger values of I^t and Ht-, given the 
relatively large mass splittings between m^' and nix- As a result, additional cuts on I^t and 
Ht (as well as p^^ and Xjj) can effectively suppress the backgrounds while keeping most of 
the signal intact, thereby optimizing the signal significance. 

Similarly, p-l^ and Meg distributions at the 7 TeV LHC after precuts are shown in Fig. [2] 
for three benchmark points {mB',mx) = (300,1), (500,1), (800,1) GeV, as well as the 
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FIG. 1: IpT and Ht distributions at the Tevatron for SM backgrounds and three signal benchmark 
points {mB',mx) = (200,1), (300,1), (400,1) GeV, after precuts. 
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Pj. and Afcff distributions at 7 TeV LHC for SM backgrounds and three signal benchmark 
mB',mx) = (300,1), (500,1), (800,1) GeV, after precuts. 



dominant SM backgrounds. The contributions from Ziyv)]], W{(iv)jj and ti are similar 
after precuts. However, for tt, the differential cross section distributions drop quickly with 
increasing ^y, pj. and Mefj , since the ti distributions are enhanced in the region below the 
t mass (for IpT and p^) ^^d ^m^ (for Mgff). With additional cuts on ^y, pj, and Mgff, the 
ti background is almost negligible. The Z{i>i')jj background, on the other hand, becomes 
dominant once additional cuts are imposed. For the signal benchmark point {mB'-,mx) = 
(300, 1) GeV, the distributions drop quickly above the mass scale of the B' . To optimize 
the cuts for such low m^/, usually no additional I^t or Pr cuts are needed and the M^s 
cut becomes the most effective in selecting the signal. For larger masses, the p!^ and Ifir 
cuts become very effective in suppressing the backgrounds. The cross sections for signal 
benchmark points and SM backgrounds after various stages of cuts are presented in the 
Appendix. 

We now determine the discovery and exclusion reach for B' at the Tevatron and the 7 TeV 
LHC. For each parameter point {mB','mx), we use the optimum cut (after precuts) that gives 
the best signal significance, with the additional requirements that S/B > 0.1 and more than 
two signal events are observed. Given the small number of signal and background events 
after cuts, we have used Poisson statistics, rather than assuming Gaussian distributions, for 



10 



Exclusion for B' B' ^ b X b X at Tevatron 



> 

0) 
O 

^00- 



250 



Discovery for B' B' ^ b X b X at Tevatron 



1 1 1 1 1 1 1 




A^\ 






"Ib = mx + 


-x^ T\\ 



200- 




150 



^~350 

> 

(U 

O 

^00- 



200- 



150- 



100^ 



50- 



- 








III 








- 


- 














^ITlg, = ITIjj + 


> 




^ 


\v 


20 fb"" 


- 




V Nidifb 


■" 


y/ 










/ 












/ 




- 




, , , , 






1 1 / /^ 





TOO 150 200 250 300 350 400 450 500 

mB,(GeV) 



ft)0 150 200 250 300 350 400 450 500 

mB.(GeV) 



FIG. 3: 95% CL Tevatron exclusion (left plot) and 3 a discovery (right plot) reach in the {mB',mx) 
plane for integrated luminosities 5, 10, and 20 fb~^. For each point in parameter space, the cut 
with the best significance has been chosen. 

both signal and backgrounds. 

Figure [3] shows the 95% CL Tevatron exclusion and 3a (Gaussian equivalent^) discovery 
contours in the {mB>,mx) plane. For relatively small values of mx, ttlb' could be excluded 
up to 440, 460 and 480 GeV, or could be discovered at 3a up to 400, 420 and 440 GeV 
for integrated luminosities of 5, 10, and 20 fb~^, respectively. For small mass splittings 
rriB' —nixi the 6-jets become soft and the amount of transverse missing energy gets smaller. 
It is more challenging to select signals out of the SM backgrounds with such soft decay 
products. This explains the gap between the exclusion/discovery contours and the dashed 
line, which corresponds to the threshold for the on-shell decay of B' — >■ hX. With 20 fb~^ 
integrated luminosity, masses nix as large as 290 GeV may be excluded, and masses nix as 
large as 270 GeV may be discovered. 

Figure m shows the 95% CL exclusion and 3o" (Gaussian equivalent) discovery contours for 
a 7 TeV early LHC run, for integrated luminosities 0.1, 1, and 10 fb~^. With just 0.1 fb~^, 
the LHC exclusion reach for m^' of 480 GeV exceeds the Tevatron exclusion reach with 
20 fb~^ luminosity. With 1 fb^^ and nix ~ 0, all regions of rra^' in the perturbative Yukawa 
coupling region can be covered. Exclusions of m^/ up about 800 GeV could be achieved 
with 10 fb~^ integrated luminosity. Note that at the LHC, we could tolerate a much smaller 
niB' — nix almost up to the on-shell decay threshold for small m^'. The 3cr discovery reach 
for iriB' is about 380, 540, 700 GeV for integrated luminosities of 0.1, 1, and 10 fb~^. The 
reach in nix is greatly enhanced at the LHC as well. It could be excluded up to 330 and 
410 GeV, or to be discovered up to 260 and 360 GeV with 1 and 10 fb~^ data. 

Note that the exclusion curve for the Tevatron in Fig. |3] fails to reach the mx = niB' — nib 
line for any values of m^', because of the small missing energy in this region. As evident in 



^ By Gaussian equivalent, we mean that we have converted the one-sided Poisson probability into the 
equivalent a deviation in a two-sided Gaussian distribution, which is more commonly used in the literature. 
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FIG. 4: 95% CL LHC 7 TeV exclusion (left plot) and 3 a discovery (right plot) reach in the 
(mB/,mx) plane for integrated luminosities 0.1, 1, and 10 fb~^. For each point in parameter 
space, the cut with the best significance has been chosen. 
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FIG. 5: Model-independent 95% exclusion, 3a and 5a discovery cross section reaches for light mx 
{mx = 1 GeV) at the Tevatron (left plot) and 7 TeV LHC (right plot). Also shown are the QCD 
pair production cross sections for B'B' (solid red curves). 

Fig. m however, this is not true at the LHC, where the energy and cross section for low-mass 
B's are so large that recoil of the B' B' system against initial state radiation jets gives events 
with sufficient IpT to pass the cuts, even for mx = rns' — Tn^,. For further discussion, see 
ReL [45i]. 

To present our exclusion and discovery reaches in a more model-independent way, in 
Fig. O we show the collider reaches of the hXhX production cross section [or equivalently, 
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a{B'B') X B{B' — t- hXY\ as a function of mB' for 95% CL exclusion, 3o" and 5o" discovery 
for various luminosities at the Tevatron and the 7 TeV LHC, with nix fixed to 1 GeV. At 
the Tevatron with 20 fb~^, production cross sections of 5—200 fb could be excluded at 95% 
CL for the mass of B' in the range of 100—500 GeV. The limits get better for higher masses 
due to the more energetic final state particles and large IpT and M^e in the signal process. 
For 5(7 discovery, the cross section reach is about 20—400 fb. At the 7 TeV LHC, with 1 fb~^ 
luminosity we could reach an exclusion limit of about 25 fb for m^' around 1 TeV. A 5(j 
reach of 20 fb can be achieved with 10 fb~^ luminosity. 

For the purpose of illustration, we also show the QCD pair production cross sections of 
B'B' (solid curves) in Fig. [51 For B'B' pair production in other new physics models, one 
can easily read out the collider reach of mB' by comparing a{B'B') x B{B' — )■ bX)"^ in those 
models with the cross section reach curves in Fig. [51 

Note that we have taken rrix = 1 GeV when presenting the cross section reaches at 
colliders. However, as evident from Figs. [31 and [H the reach in ttlb' is almost independent of 
rux for small and moderate values of mx, unless the mass splitting of m^' ~ ^^x becomes 
small. 



VI. CONCLUSIONS 

In this study, we have considered the possibility of pair production of new charge — | 
quarks B' that decay directly to 6-quarks and long-lived neutral particles X. The resulting 
signal is B'B' — )■ hhlfij'-, which is common to many new physics theories, as discussed in 
Sec. [Til Since the bblpT signal is common to many models, we have also presented detection 
prospects in terms of the pair production cross-section a{bbXX), for varying m^/ and various 
integrated luminosities. This analysis thus accommodates many models in which there are 
new contributions to B' B' production, as well as models where B' is not spin-|, as in the 
case of bottom squarks.'^ 

We have estimated the sensitivity of the Tevatron and LHC to new physics resulting in 
this signal. Currently published bounds have been summarized in detail in Sec. IIIIj very 
roughly, however, and translating bounds on related processes to the case of B'B' production, 
the current limits are m^/ ^ 370 GeV. From our analysis, we expect that these results may 
be improved to uib' ^ 440 GeV for small mx-, given our optimized cuts and 5 fb^^ of 
data, which is currently available at the Tevatron. We also find that additional Tevatron 
data will marginally improve this bound: for 20 fb~^ of data, models with uibi < 480 GeV 
{rnx ^ 150 GeV) can be excluded at 95% CL. 

These results may be further improved at the LHC with an analysis of 1 fb^^ of data, 
which has already been accumulated. With 1 fb~^ of integrated luminosity, the LHC physics 
run can probe any models with m^/ ^ 620 GeV, provided mx ^ 250 GeV. Models with B' 
quarks that get mass from electroweak symmetry breaking are bounded by the requirement 
of perturbative Yukawa couplings to have masses below this mass. Early LHC data may 
therefore probe the full range of possible quark masses in these models. In particular, the 
early LHC will probe all WIMPless models that could explain the data of DAMA and 
CoGeNT (assuming dominant coupling to 3rd generation SM quarks). Even null results 



If B' is not spin-i, then the angular distributions of the outgoing jets wih change, altering cut efficiencies. 
In this case, the comparison is only approximate. 
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from the search discussed here will therefore be of significant interest. 

These same detection prospects are, of course, also applicable to little Higgs models, or 
asymmetric dark matter models arising from hidden sector baryogenesis. For these models, 
however, theoretical considerations provide no expected upper bound on rriB'- The mass 
reach which the LHC can achieve with greater luminosity is thus of interest. With 10 fb~^ 
of data, the LHC can probe asymmetric dark matter models and other similar frameworks 
at the 95% CL for ttib' < 800 GeV, provided mx < 200 GeV. 

UED models have perhaps the least constrained theoretical motivation, since the B' 
mass is not bounded by Yukawa coupling perturbativity and relatively small B' — X mass 
splittings are perfectly plausible. With 10 fb~^ the LHC can probe models with uix as large 
as 410 GeV at the 95% CL. This maximum reach is obtained for uib' in the 600 — 700 GeV 
range. 

All of these detection prospects can be easily translated into mass reaches for bottom 
squarks decaying directly to hx\- With 1 fb~^, the LHC can probe models with m^ < 
400 GeV at 95% CL (provided m^o < 150 GeV). With 10 fb-\ the LHC reach increases to 
mi < 520 GeV. 

It is worthwhile to compare the mass reach of this B' search to that of the T' search 
examined in Ref. [l|. There it was found that, with 1 fb~^ of data, the LHC could probe 
low-mass dark matter models at 2)a for m^' < 490 GeV. The B' search described here 
has similar reach for the same luminosity and required signal significance. However, the 
T' detection prospects were seen to drop rapidly with increasing mx, with no sensitivity 
at all expected for mx > 180 GeV. In contrast, the detection prospects for this B' search 
are almost unchanged for mx ^ 200 GeV (assuming 3cr significance). This difference is 
attributable to the large mass of the top quark; for relatively heavy X, there is very little 
phase space left for the T' — )■ tX decay. Although the tt^r signals provides many more 
handles, in the end, the naive expectation holds true: the reaches in ms' and m^T' are roughly 
similar, and for a fixed new quark mass, the dark matter mass reach of the B' search exceeds 
that of the T' search by roughly m,t — m,h ^ 170 GeV. 

The analysis presented here determines the prospects for detecting an excess in events 
with fe-jets and missing Et- It is more difficult to determine if the excess arises from the pair 
production of B', decaying via B' — )■ bX . To determine the masses m^/ and mx would be 
harder still. It would be interesting to determine the prospects for the LHC to make these 
measurements. 
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Appendix: Impact of Cuts on Signal and Backgrounds 

In this Appendix, we present tables listing the cross sections after cuts for the B'B' signal 
and the main SM backgrounds (Tables HllITl). In the upper section of each table, each line 
gives the cross section after including all cuts above. In the lower section, each line gives 
the cross section after including the cut on that line, and all precuts. For the signal, three 
examples with nix = I GeV and m^/ = 200, 300, and 400 (300, 500 and 800) GeV were 
chosen for the Tevatron (7 TeV LHC). The W and Z cross sections in parentheses were 
simulated with a cut on I^t > 20 (60) GeV for the Tevatron (LHC) and at least 2 jets in 
the parton-level generation. 



TABLE L Signal and background cross sections in pb after cuts for signal and dominant back- 
grounds at the Tevatron. The signal examples are for mx = 1 GeV and niB' = 200, 300, and 400 
GeV as indicated. The W cross sections in parentheses were simulated with a cut on ^t > 20 GeV 
and at least 2 jets in the parton-level generation. From the 160 independent combinations of final 
cuts used for the cut optimization, the three cuts that optimize the significance for these three 
mass points are displayed in the table. Momenta and masses are in GeV. 





B' (200) 


B' (300) 


B' (400) 


W^+jets 


Z — )• z^z^-|-jets 


ti-|-jets 


No cut 


2.62 


0.195 


0.0154 


(632.45) 


(21.103) 


5.628 


leptons 


2.24 


0.169 


0.0134 


(229.22) 


(16.516) 


2.365 


2 < jets < 3 


1.89 


0.143 


0.0109 


(33.80) 


(7.962) 


0.456 


aj,j^ < 165° 


1.66 


0.125 


0.0097 


(29.35) 


(7.171) 


0.362 


^T >40 


1.55 


0.122 


0.0096 


17.05 


5.221 


0.235 


^^ > 80 - 40 X A,^S^'^^*^) 


1.52 


0.121 


0.0095 


16.45 


5.042 


0.207 


A(/>(^r,jets) > 0.6 


1.43 


0.112 


0.0086 


16.29 


4.957 


0.188 


^=SS^-t 


1.42 


0.111 


0.0086 


15.89 


4.869 


0.179 


X,j^{p'^+P>^)/Ht> 0.75 


1.25 


0.102 


0.0081 


13.48 


4.197 


0.079 


p^^ > 20 


1.25 


0.102 


0.0081 


13.48 


4.197 


0.079 


Ht >60 


1.25 


0.102 


0.0081 


10.56 


3.680 


0.078 


> 2 6-jets, 6-jet hardest jet 


0.43 


0.035 


0.0026 


0.11 


0.037 


0.018 


All precuts 


0.43 


0.035 


0.0026 


0.11 


0.037 


0.018 


Xjj > 0.9, ^T > 40, Ht > 300 


0.018 






8.59 • 10-5 


1.11-10-^ 


1.68 • 10-^ 


Xjj > 0.9, ^T > 150, Ht > 300 




0.0043 




4.56 • 10-5 


7.40 • 10-5 


4.11-10-5 


Xjj > 0.9, ^T > 250, Ht > 300 






5.30 • 10"^ 


1.95 • 10-5 


4.02 • 10-5 


1.39 • 10-5 
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